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Abstract 

We consider mass-conserving self-similar solutions for Smoluchowski's 
coagulation equation with kernel K(^,rj) — (£,r]) x with A £ (0,1/2). It 
is known that such self-similar solutions g(x) satisfy that x~ 1+2X g(x) is 
bounded above and below as x — > 0. In this paper we describe in detail via 
formal asymptotics the qualitative behavior of a suitably rescaled function 
h(x) = h\x~ 1+2X g(x) in the limit A — > 0. It turns out that h ~ 1 + 
Cx ' cos (VAlog x) as x — > 0. As x becomes larger h develops peaks of 
height 1/A that are separated by large regions where h is small. Finally, h 
converges to zero exponentially fast as x — > oo. Our analysis is based on 
different approximations of a nonlocal operator, that reduces the original 
equation in certain regimes to a system of ODE. 



1 Introduction 



In this article we consider self-similar solutions to Smoluchowski's mean-field 
model for coagulation [H [TB] that is given by the equation 

1 ft 

dtf(£,t) = - / K(£-ri,Ti)f(£-Ti,t)f(ri,t)dri 

, (1-1) 
K(t,fl)f{t,t)f(v,t)dri=:Q\f\ , 







where f(£,t) denotes the number density of clusters of size £ at time t. The 
kernel K(£,r/) describes the rate of coalescence of clusters of size £ and rj and 
subsumes all the microscopic aspects of the coagulation process. It is well- 
known, cf. for example the review article [TU] and references therein, that if the 
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kernel grows at most linearly, then the initial value problem for (jl.ip is well- 
posed for initial data with finite mass and the mass is conserved for all times. For 
homogeneous kernels it is furthermore expected that solutions converge to a self- 
similar form for large times. In fact, for the kernels K(£, rj) = 1 and K(£, ij) — 
£ + r], the constant and additive kernel respectively, this issue has by now been 
completely solved. It has been established [U [T3] that besides explicitly known 
exponentially decaying self-similar solutions also solutions with algebraic decay 
exist. In [T3] their domains of attraction could also be completely characterized. 
However, not much is known about self-similar solutions for other kernels than 
these solvable ones. Existence of fast-decaying self-similar solutions has been 
established for a large class of kernels [S1[S], and local properties of such solutions 
have been investigated in [3j [6l [9l [14], but it is not known whether they are 
unique in the class of solutions with finite mass. It is also not clear, not even 
on the formal level, whether self-similar solutions with algebraic decay exist. 

In this article we will consider mass-conserving self-similar solutions to coag- 
ulation equations with the following product kernel of homogeneity 2A € (0, 1): 

k(h,v) = ({;v) x , °< A <^ (1-2) 

Existence of a solution for this kernel has been established in [8]. The purpose 
of this paper is to give an asymptotic description of such a solution in the regime 
A ->■ 0. 

Let us describe briefly what is known for self-similar solutions for kernels as 
in (|1.2j) and compare it to results for the so-called sum kernels. To fix ideas, we 
restrict ourselves to kernels of the form 

V) = fV + £V , " + P = 2A, < a < 0, (1.3) 

even though most of the results that we mention also apply to more general 
kernels with the same growth behavior as the ones in (| X . 3|) . We will need in 
particular to distinguish between the case a — 0, the sum kernel, and the case 
a > 0, the product kernel. It has long been predicted [HJ [T7] that self-similar 
solutions for kernels as in (|1.3[) exhibit a singular power-law behavior of the 
form x~ T with r < 1 + 2A in the case a — 0, and r = 1 + 2A for the case a > 0. 
This has been rigorously proved for the case a = in [9] and for the case a > 
in [6] [14]. As has been pointed out in [9], the next order behavior for small 
clusters in the case a — can then easily be established and is as predicted 
by the physicists. In the case a > 0, however, the next order behavior has not 
been known [17] , and only numerical simulations suggested that it is oscillatory 

[am- 
our goal in this paper is to describe formally how to construct mass- 
conserving self-similar solutions for kernels as in (ll.2j) in the limit A — > and 
by this also describe their asymptotic behavior on the whole positive real line. 
While we believe that such self-similar solutions are unique, our approach only 
gives the construction of one such solution. In forthcoming work |15] we will 
also show how to make this construction rigorous. We will see that in the limit 
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A — » the oscillatory character of the solutions becomes very explicit and can be 
interpreted in terms of simple ODEs coupled with explicitly solvable equations. 

In the limit A — > + the product kernel becomes close to the sum kernel 
x 2X + y 2X . As described above the power law for the sum kernel is different from 
the one for the product kernel. This is due to the fact that, thinking in terms 
of coagulating particles, the physical behaviour of these particles is different in 
these two cases. Indeed, the power law obtained for the sum case shows that 
particles with small x coalesce mostly with particles that are much bigger than 
themselves. On the contrary, in the case of the product kernel small particles 
interact mostly with the ones having a comparable size. The point A = can 
be thought as a bifurcation point where both kind of behaviours take place 
simultaneously. 



2 Preliminaries and Overview 
2.1 Equation for self-similar solutions 

We are now going to derive the equations that are solved by self-similar solutions. 
It is known that the coagulation equation 1)1.10 can be written in the follow- 
ing conservative form that makes the conservation of the number of monomers 
/o°° %f (£' d£, transparent. 

dt (£/ & t)) + dt (^f* dr, jH d pK fa, p) v f ( V , t) f (p, t^j = . (2.1) 
Using the self-similar variables 

f&t)=t-*$xg(x), x = -^, (2.2) 

t 1-2A 

as well as the form of the kernel (|1.2[) . equation (|2.1[) becomes 

~d x (j^xj) + d * d V dzK (y, z) yg (y) g (z)^j = . 

Integrating this equation, assuming decay of the solutions as x — > oo and im- 
posing absence of particle fluxes, we obtain 



dy dzK (y,z)yg(y)g(z) . (2.3) 



x 2 g (x) 
(1-2A) ■ 

The balance of the terms on both sides of (|2.3[) suggests the following power law 
behaviour for g near the origin, 

g(a) ~ (l-2A)x^ as ^ 0+ < H >= B{l-il-X) {2A) 
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where B (•, •) is the classical Beta function PQ. For further reference we notice 
that 

B (1— A, 1— A) ~ 1 + 2A + o(A) asA^O. (2.5) 
In order to remove the power law behaviour we introduce the new function 

h ( x ) = -—jr~^ xl+2X 9 ( x ) • (2-6) 

Then ft solves 



h{x)=H x x zx - 1 dyy- 2X h{y) dzK (y, z) z - {l+lx) ft (z) . (2.7) 

Jo Jx-y 

In the rest of the paper we will study solutions of (|2.7[) in the limit A — > 
that satisfy — ► as a; — > oo. Notice that (12.7[) has the explicit constant 
solution ft, (x) = 1 that corresponds to the well-known solution with infinite 
mass g (x) = nzffxi ^1+^ m tne formulation ()2.3|) . 

2.2 Reformulation as a Volterra integro-differential equa- 
tion 

A main idea in our approach is to reformulate the problem of finding solutions of 
(|2.7[) as a shooting problem for a Volterra integro-differential equation. Indeed, 
we can rewrite (|2.7I) as 

h (*) = T77^TT U (») « (*) + ^a^ 1 f f dzy" A z- (A+1) ft (y) ft (z) , 

(2.8) 

u(x) = {l-\)x x - 1 dyy- x h(y) , v (x) = Xx x dzz'^+^h (z) . (2.9) 
Jo ii 

Differentiating (|2.9[) we obtain the equations 

ra x = — (1 — A) u + (1 — A) h , xv x = Xv — Aft . (2-10) 

Thus the solution of (|2~7) solves (|2~T0|) with ft as in (12.81) . that has the advantage 
of being a Volterra integro-differential equation. However, the problem (|2.8p . 
(|2.10| admits solutions that in general do not decay as x — > oo and therefore 
they do not provide solutions of (|2.7|) . Therefore, among all the solutions of 
|, (|2.10p we must select those satisfying 

v(x)~^Q, ft(x)->0, u(x)^0 as x ->• oo . (2-11) 
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2.3 Overview of different regimes 



We now sketch the asymptotics of the solution that we are going to construct. 
It passes, roughly speaking, through three stages that are described in detail in 
Sections O 2] and [5] respectively. 

Section [3] discusses the behavior of the solution near x = 0. As pointed out 
before, near the origin the solution is oscillatory and behaves as 

cos (a (A) log (x) + ip) as x — > + (2-12) 

where /3(A) ~ A/2 and a(A) ~ ^/X as A — > 0. Our analysis consists of con- 
structing a solution by starting with h as in (j2 . 12[) and using K — Ce ttp as a 
shooting parameter. However, since equation (|2.7p is invariant under the rescal- 
ing x n- ax for a > 0, we can restrict the range of K to [1, exp 2n((3/a)) (see also 
the comment in Section [3.2[1 . In this regime near the origin we can approximate 
(|2.8I) and (|2.9p by a nonlinear ODE (see (|3.19[1 ). that is a perturbation of a 
simple ODE system. For the perturbed system we can use an adiabatic approx- 
imation to compute the increase of an associated energy E along trajectories. 
As we discuss in Section 13.51 this first ODE approximation is valid as long as 
E <C 1/A. When E ~ 1/A we enter a new regime that is described in Section UJ 
In this regime h develops peaks of height 1/A and width of order one that are 
separated by wide regions in which h is small. The regime where h is small is 
again described by an ODE (cf. Section |4~2]) . while the peaks are described by 
an integro-differential equation (see Section l4~Tj) . The analysis of these respec- 
tive regimes is done in Sections 14.31 and 14. 4[ while their coupling is described in 
Section [4.51 In Section [3] we will then show by a continuity argument that there 
exists a shooting parameter such that the corresponding solution h converges 
to zero as x — > oo. In the third regime, that is discussed in Section [6l h then 
decays exponentially fast to zero. 



3 The behaviour as x — > 
3.1 Oscillations 

It has been observed in [?[ [TT] that the self-similar solutions associated with 
(|1.1[) exhibit oscillations for small values of x. This oscillatory behaviour can be 
seen as follows. If we make the ansatz h(x) = 1 + ax^ + ■ ■ ■ with some /i G C, 
plug this into (|2.7j) and keep only the leading order terms, we obtain that /x 
must satisfy 

i= tdyy-^ r dzK(y,z)z-<- 1 +* x >+» 
"A Jo Ji-y 

+ / dyy- 2X +» / dzK(y,z)z^ 1+2X K 
Jo Ji-y 

After some computations, we find that this is equivalent to 

* A (^) := (A-/u) B (1-A, 1-A) - (2A-^i) B (1-A, l-A+£t) = . (3.1) 
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One can show that ty\ has two complex conjugate roots fj,f — ±a (A) i+fi (A) 
with positive real part and nonzero imaginary part. Indeed, in the limit A — > 
this can be seen easily as the dominating terms in (I3.1[) are, assuming 3> A, 

a (n) ~ (A-/i)(l+2A) - (2A-/x)(l-/x+2A + /i 2 ) 

and thus the roots of ^ \ satisfy 

fif = ~±VXi + o(x§^) as A^0+. (3.2) 

Consequently we expect the following asymptotics for the solutions of (12. 7[) : 

cos (a (A) log (x) + ip) as x — > + (3-3) 

for suitable real constants C and <£>. 

Notice that the asymptotics (|3.3[) indicate that two consecutive local maxima 

of h, called x\ and X2, satisfy ~ e xp f° r small A. In comparison, 

to double the amplitude of (|3.3I) we must multiply i by a number of order 

exp ^ 21 °s( 2 ) ^ . Thus, the behaviour Q3.3P is essentially oscillatory, with a growth 

in the amplitude that takes place on a much larger scale. 

3.2 A new set of variables 

For our forthcoming analysis it will be more convenient to use the following 
variables: 

x = e x , y = e Y , z = e z , (3.4) 
h(x)=H(X), u(x) = U(X), v(x)=V(X). (3.5) 

Then (1211), (|2T0| becomes 

H= J ^-UV + I[H] , (3.6) 
^ = -(l-X)U+(l-X)H, ^ = XV-XH, (3.7) 

where 

I [H] (X) =H X [ dYe {1 - x ^ Y - x ^H (Y) f dZe~ x ( z ~ x>> H (Z) . 

J — oo J\og(e x —e Y ) 

(3.8) 

The system (|3.6[) - (|3.8p has two explicit solutions, namely (H,U,V) 
= (1,1,1) and (H,U,V) = (0,0,0). Our goal is to construct a solution of 
(I3.6p - (|3.8I) for small values of A with the property 

lim (H (X) ,U (X) ,V (X)) = (1,1,1) , (3.9) 

X^t — oo 

lim (H (X) , U (X) , V (X)) = (0, 0, 0) . (3.10) 
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We will see that a solution that satisfies (|3.10p decays exponentially fast to zero. 
This in particular implies that the corresponding self-similar solution g has finite 
mass. 

Equations (|3.4j) . (|3.5[) and (13.31) yield the following asymptotics of the solu- 
tions as X — > — oo: 



(H (X) , U (X) , V (X)) ~ (1, 1, 1) + Re [K 1, 



(1-A) 



-A 



(1-A)+ M +' (/i+-A), 

(3.11) 

where /x+ is as in and if = Ce^ e C. 

The complex number K determines the solution of (|3.6I) - (|3.8[) uniquely. How- 
ever, we notice that (I3.6l) - (|3.8p is invariant under translations X — > X + a and 

thus all the complex numbers in the spiral S\ t x — |ife M+a : a € C C yield 
the same solution up to translations. Thus we can identify the solutions of 
(|3-6[) - (|3.8[) satisfying (|3.9[) . up to translations, with the set of positive real num- 
bers contained between two consecutive intersections of S\ with the real axis. 
In particular, given K = 1 we have, since ^ = a(X)i + (3 (A), that the next 

consecutive point in IR + C\Sx,i is exp ^ 2 a(A)^ ) ■ Therefore, there is a one-to-one 

and the 



correspondence between the points in the interval 

solutions of dSU-dSU) satisfying (j3~9| . 

The main result of this paper is that we show, using asymptotic arguments, 



that for small A a value of K € 



exp ( 2 ^a(x^ ) ) exists such that (|3. 10[) holds. 




Figure 3.1: Oscillatory behavior of U (left) and phase plane for U, V (right). 



3.3 The ODE regime 

The key idea in computing the asymptotics of the solutions of (|3.6|) - (|3.8[) is to 
obtain suitable approximations of the operator / [H] (X) in (|3.8I) for small A. 
The formula for H\ (cf. (|2.4p ). combined with (|2.5j) and (|3.8j) . suggests the 
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approximation 

/x 
e {Y - x) H (Y)dY / H{Z)dZ. (3.12) 

-oo J\og{e x — e Y ) 

We will discuss the consistency of this assumption in Section [7] 

The operator I [H] (X) can be further approximated if Y/X changes signifi- 
cantly faster than H(Y)/ H(X), as we saw is the case for X — > — oo. Then it is 
natural to approximate H{Y) by H(X) and we obtain 



I[H](X)~ \H(X) f e {Y ~ x) dY f H(Z)dZ 

J-oo Jlofr(e x -e Y ) 



and changing the order of the integrals 



f x 

I [H] (X) ~ XH (X) / H(Z)e z - x dZ. (3.13) 

J — OO 

From (|3.7|) we have 

/X pX 
H{Z)e {1 - x){z - X) dZ ~ / H{Z)e {z - x) dZ 
-oo J —oo 

which together with (|3.13p leads to 

I[H]~\HU. (3.14) 

Using (|3.14l) as well as (|2.4[) and (I2.5[) we obtain the following approximation of 
(|3.6[) up to order A: 

H = UV + XUV (U - 1) (3.15) 
and plugging this approximation into (|3.7|) , p.8[) we obtain 

^ = - (1-A) U + (1-A) C/y + A (1-A) C/y (17—1) , (3.16) 

— = XV - XUV - X 2 UV{U - 1) . (3.17) 
dX 

As (|3.2[) and (|3.1ip imply that the changes of (U— 1) compared to those of (U—l) 
are of order -\/A, we make the change of variables 

y = i + VAw , £ = Va~a:, (3.1s) 

that transform (|3.16j) . (|3.17|) up to order a/A into 

^ = Uu + VXU (U - 1) , ^ = 1 - C/ + VAw(1 - {/) . (3.19) 
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3.4 Analysis of the ODE (13A91) 

Let us denote by (U\, oj\) any solution of (I3.19[) . We first notice that the asymp- 
totics of the solutions of (|3.19p as (U\,L)\) — > (1,0) agree with those obtained 
in p. lip . Indeed, in the limit A — > 0, we can compute the asymptotics of (U, uj) 
as f -» -oo using ([57TT]) . (|3~2]) and p~T5)) : 

((C/ A -l)(0,^(0)^Re(^(l, l )e( l+ °(^))«) as£->-oo. (3.20) 

Linearizing (|3.19p near (U,u>) = (1,0) we obtain the same asymptotics for the 
functions (U\ (£) , ui\ (£)) as £ — > — oo. This gives the desired matching between 
the solutions of the linearization of (|3.6p ~ (13.81) around (if, £/, V) = (1,1,1) and 
the solutions of the approximated problem (|3.19p . 

The approximation (|3 . 20[) is valid as long as \(U\ — 1 (£) ,ui\ (£))| is small. 
However, a detailed analysis of the nonlinear problem (|3.19l) is needed if 
\(U\ — 1 (£) ,lu\ (£))| becomes of order one. In order to describe the solutions 
of (|3.19p in this regime, we notice that this equation is a perturbation of the 
equation 

f = ^». ^ 

This equation can be explicitly integrated, since the following quantity is con- 
served along trajectories: 

E = - log (U Q ) + (U Q - 1) + ^ . (3.22) 

The solutions of (|3.2ip are periodic, but this behavior is not compatible with 
the exponential growth obtained in (|3.20p . This growth is due to the increase in 
E produced by the terms of order VA in (|3.19p . We now compute this change 
of energy for values of (U, lj) of order one. Using Q3.19D and (|3.22p we obtain 



(Ux-lf +u 2 x (l-Ux) . (3.23) 



Since the solutions of p,19p are close to those of p.2ip during finite time 
intervals, we can adiabatically compute the change in E. More precisely, the 
solutions of (I3.2ip satisfying p.22p are periodic with a period T (E) given by 

T (E) = y/2 / 1 (3.24) 

Ju-(E,o) r]y/E + \og(rj) - (77- 1) 

where the functions U- (E, u) < U + (E, uj) are defined as the roots of the equa- 
tion 

Iog(C0-(^-l) = y-^ (3-25) 

for any given value of uj and E > ^- . Integrating (I3.23P and using the adiabatic 
approximation, we can then approximate the change of E during an interval of 
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length T (E) as 

e(t { e)j_e(o) = jnm ^ ^ _ i)2 + ^ ^ (i _ ^ ^ ^ = $ ^ 

(3.26) 

where (C/o (£) , wo (£)) is a solution of (|3.2ip satisfying p. 221) . 

It turns out that the function $ (E) is positive for any E > 0. This follows, 
using the second equation in (|3 . 2 1 1) . via 

/\/2E rV2E 
_ [(1 - E/_ (E,lj)) +uj 2 ] doj - / _ [(1 - U+ (E,uj)) + lj 2 ] du, 
-V2E J-V2E 

whence 

$ (E) = 2 / (U+ (E, u) - U- (E, Lo))duj>0. (3.27) 







We can compute the asymptotics of (|3.27l) as E — > oo using (|3.25l) . No- 
tice that the leading contribution to the integral is given by U+ [E, u>) , since 
U- {E,uj) < 1. Then 

<$>{E)~2j \E - — \ duj = -7j-E§ as £^oo. (3.28) 
3.5 Range of validity of the ODE regime 

We have obtained that, as long as the approximation (|3.19l) is valid, we can ap- 
proximate the functions (U\ (£) , oj\ (£)) by (/To (^),uo (£)) while E is computed 
via the iterative recursion (|3 . 26[) . Due to (|3.28p these values of E increase an 

amount of order vA(l + E) 1 in each period of size T (E) . We now proceed to 
discuss the range of validity of the different approximations that have been used 
to derive and QMB- 

To examine the range of values of E for which the adiabatic approximation 
is valid, we need to show that the value of E (£) remains close to E (0) for 
< £ < T (E) . Due to (|3.25l) and (|3 . 2T[) the adiabatic approximation is valid 
as long as 

sup \E(£) -£(0)| < E(0) . (3.29) 

0<{<T(_B) 

Let us remark that (|3 . 29[) does not follow immediately from the inequality 
\E (T (E (0))) - E (0)| < E (0) because the right-hand side of (QEH) does not 
have a sign and therefore it is not clear that 

sup \E(O~E(0)\<C\E(T(E(0)))-E(0)\ (3.30) 

0<i<T(E) 

for some C > 0. A priori it is not possible to rule out the possibility of big 
changes of E (£) along each cycle balanced in such a way that the overall change 
of E along the cycle is just $ (E) . However, as shown in Appendix [S] it turns 
out that (|3.30[) holds true. 
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Using (I3.26[) , (13.28)) and (|3.30[) it follows that the adiabatic approximation 
is valid as long as vA(l + E) 1 <C (1 + E) and this is satisfied as long as 



A (! + £)< 1. 



(3.31) 



We now study the range of validity of the approximation (|3.14|) that is the 
main ingredient in deriving the approximate problem (|3.19|) . The main assump- 
tion used in the derivation of (13. 14)) is that the characteristic length scale for 
which (H (X) — 1) has changes comparable to itself is much larger than one. 
Due to (|3.15[) it follows that as long as V remains close to 1 the characteristic 
length scale for (H — 1) is the same as for ([/ — 1) . Notice that (V — 1) is small 
if y/\ui <C 1. Since we have < E, the condition VAa; <C 1 holds if f|3.31[) 
is true, and so the characteristic length scales for {H — 1) and (U — 1) are the 
same under the assumption (|3.31l) . For \(U, uj)\ of order one, the evolution of 
(|3 . 2 1 1) takes place in the length scale £, or equivalently for changes in X of order 

^> 1. Therefore, the condition in (|3.31[) that allows us to obtain (|3.14|) is 
immediately satisfied and we can restrict our analysis to the case \(U, uj)\ ^> 1. 
The first equation in (|3.21[) combined with (|3.18p implies that 



The left-hand side measures the relative variations in U with respect to 
changes in X. The characteristic length scale that describes the changes in this 
quantity is large as long as \/A <C 1 which is again true if p.31[) holds. 

Therefore, if the condition p. 311) is satisfied, we can apply simultaneously 
all the approximations that lead to (13.19)) , the adiabatic approximation yielding 
13.26)) as well as the condition | V — 1 j -C 1. However, the three assumptions 
break down simultaneously if E becomes of order j. Since (|3.26p and (|3.28|) 
imply that E increases to arbitrarily large values, the failure of (|3.31[) happens 
for every solution of (I3.6[) - (|3.8[) satisfying (|3.11[) . 

4 The intermediate regime 

In order to describe the solutions of p.6[) - (|3.8|) when XE w 1 we introduce a 
new group of variables. It turns out to be convenient to split the regime into 
the one where U (or H) is of order one or smaller, and the one where U (or 
H) is large. Notice that (|3.25j) implies that for a given value of E ~ \ the 
minimum of U scales as e -1 / A , and the maximum scales as 1/A. We remark 
that the forthcoming analysis will show that as long as V remains of order one 
and is not small the values of H and U have the same order of magnitude. This 
assumption will be made implicitly in all the remaining computations and will 
be justified by the self-consistency of the derived asymptotics. 



dX 



(3.32) 
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4.1 The integro-differential equation regime 

We begin by studying the case in which U (and H) are large. Since we are 
interested in the case E rs ~ , equation (|3.25[) suggests the rescaling U ~ j , w « 
-j=. Then (|3.32[) yields a characteristic length scale for X of order one. This 
suggests introducing the new set of variables 

U{X) = ±U(X), V(X)=V(X), H(X) = jU(X). (4.1) 

Plugging dUI]) into (|3"1) ]| -([3~% ]I we obtain 

H = ^UV + ^I[n}, (4.2) 

^ = -(l_A)W+(l-A)K, ^| = AV -^' ( 4 - 3 ) 

where the operator / [H] is as in ()3.8|) . Using the approximation (|3.12j) that will 
be seen to be still valid in this region we obtain 

1 1 [%] (X) = (1 + (A)) / e^ Y - x ^n (Y) dY [ H (Z) dZ . 



A 

Taking the limit A -> in (f4~2]) -(f4~3 l) we obtain 



log( e x_ 



X_„Y\ 



e {Y - x) H(Y)dY H(Z) ,dZ (4.4) 

-oc J log(e x — eX ) 

4.2 The ODE regime 

If XE s» 1 but Z7 (and i?) remains bounded we can approximate (|3.6p ~ (I3.8[) as 
follows. Keeping just the leading terms in (|3.6|) . (|3.8|) we arrive at the approxi- 
mation H = UV. Plugging this approximation into p.7[) we obtain 



§ = -U + UV , %= H V-UV). (4.6) 
4.3 Explicit solution of (H3D, (1431) 

We are going to compute the solutions of (14.41) . (14.5[) explicitly using Laplace 
transforms. In order to bring (|4.4j) . (|4.5p into the form of a convolution equation 
we introduce the set of variables 

x = e x , % (X) = h(x) , U(X) = u(x) , V (X) = v (x) (4.7) 

that transforms (14.411. (14.51) into 



1 r - f x h(z) 

h(x)=u(x)v(x) + - h{y)dy -^-dz , (4.8) 



xu x = —u + h , xv x = — h. (4-9) 
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Notice that this system of equations can be obtained from (j2.8[) , (|2.10[) by 
means of the rescaling h — j, u — j, v = v and taking the limit A — > + . We 
can compute explicitly a family of solutions of (|4.8[) . (|4.9|) that will be used to 
describe the solutions of (|3.6p - (|3.8l) in the limit A — > + . 

Theorem 1 For any given a £ (0, 2) and k > there exists a solution to {4.8\ ), 
|^.ff| ) i/iai satifies 

h(0) — h(oo) = u(0) = u(oo) — , lim v fx) = 1 + a, lim v (x) = 1 — a . 

(4.10) 

It is given by 

v(x)^ — I 
2m J 1 

where 7 = {7 (s) : s £ M} is a contour in the complex plane contained in the 
domain {arg (£) £ (—6*0,^0)} satisfying arg(7(s)) — > — 0q as s —00 and 
arg (7 (s)) — > 9a as s —> 00 with \ < 0$ < min 7r} . 

// < a < 1 we /iaue w x (a;) < 0, h (x) > and u (x) > for any x £ R + . 

The following asymptotics hold: 

» z M~- r n i a as x^0 + , a£(0,2) (4.12) 
1 (a) x l a 

2a 1 

«xW r / as 33 ^ 00 ' a e (°' 1 ) u I 1 ' 2 ) 



kx 



(4.11) 



/i(x) 
h(x) 



2an 
T(a) 



: a as x -> 0+ , a € (0, 2) 



(4.13) 



2 a 



1 



«r (—a) x a 



as x — >• 00 



o€(0,l)U(l,2) 



2an 

u{x) x a as x^0+ , 6(0,2) (4.14) 
1 (a) (a + 1) 

* (a) ~- «r(J) O (i-<0 ^ « , ae(0,l) 

J/a = l: 

u far) = 2e~ ra , h far) = 2 (kx) e" Ka; , u fx) = — (l - (1 + kx) e~ KX ) (4.15) 

KX 

Remark 2 The solutions with different values of k are essentially the same up 
to rescaling. Notice also that \4 . 1 0| ) implies that v (x) becomes negative for x 
sufficiently large if a > 1 . 

Remark 3 Some comments about the choice of the range of values of a are in 

1 

order. We need the restriction < a < 2 to avoid the singular points (—K) a 
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crossing the imaginary axis. It would also be possible to include the value a = 1 
in the second formula of \4-12ty and 1 Sty since the right-hand side vanishes. 
However this would result in asymptotic formulas like f ~ that do not have a 
precise meaning. Finally the constraint a < 1 in the second formula of \4-ljty 
is strictly needed, since the asymptotics of u (x) for a > 1 is proportional to j. 

Remark 4 It is interesting to note that the solutions described in Theorem [JJ 
for < a < 1 are just the self-similar solutions with algebraic decay that have 
been obtained for the coagulation equation with constant kernel in '131. This 
correspondence can be seen as follows. The solutions we obtained in Theorem^ 
satisfy 

f°° _ _ 1 f x 

v(x) = (l-a)+ g(z)dz , u(x) = - yg{y)dy 



with g (x) — MiEl, Using this formula to eliminate u, v from we obtain 

f-X px />oo 

x 2 g(x) = (i-a) yg(y)dy+ yg(y)dy g(z)dz. (4.16) 



Differentiating \4-lb^ty we obtain the equation that is satisfied by the self- 
similar solutions of 

d T F = Q[F]- (1 ~ g) f (4.17) 

T 



having the form f (£, r) = ^g \JL\ with Q [•] as in il.l]) . Equation \4- ITty can be 

transformed into 1 1.1]) using the change of variables F (£, r) = (r) a -* / (£, t) 
and t = — and correspondingly the solutions obtained in Theorem [JJ are trans- 
formed into those obtained in US]/ if < a < 1. 



Proof. Integrating the first equation in (|4.9j) we arrive at 



u(x) = - I h(z)dz. (4.18) 



x 



Using the second equation in (|4.9p we obtain 

h (z 



dz = -v (x) + v (x - y) , (4.19) 



h(z)dz = — / zv z (z) — — xv (x) + / v (z) dz , (4.20) 
Jo Jo 

assuming that all the integrals involved are convergent. Therefore, using (|4. 18[) 
and (|4T20|) we find 

1 f x 

u(x) = -v(x) + - v(z)dz. (4-21) 
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Eliminating h and u from 
(14181) . dH9l) , (rail , implies 



-XVz (x) 



/ / nn2 V (x) 

(u (x)) 2 + ^ 



, the second equation in (|4.9p . as well as 

v (z) dz 



yv y (y) v(x-y)dy + 



v (x) 



y% (2/) dy. 



Integrating by parts in the last equation we obtain 



x 2 v x (x) 



yv y (y) v(x-y)dy. 



(4.22) 



Equation (|4.22j) is a convolution equation. In order to solve it we introduce 
the Laplace transform of v (x) 



(C) = / v (x) e~ Cx dx . 



(4.23) 



Then (|4.22l) becomes ^j(C w ) + W ^(C W ) — 0- This equation can be integrated 
explicitly. Its only solution that does not have singularities and takes real values 
along the line £ £ R + is 



w(C) = 



1 



C - K 
( a + K 



Inverting the Laplace transform we obtain 



C a - K 
( a + K 



(4.24) 



(4.25) 



The change of variables (x = ( transforms this formula into (|4.11j) . The 
choice of the contour of integration ensures the exponential convergence of the 
integral as well as the absence of singularities of the integrand along the contour 
of integration. Taking the limits x — > + and i->oowe obtain (|4.10[) . 

The negativity of v x for < a < 1 can be proved by differentiating (|4.25|) 
and deforming the contour of integration to a double half-line following M _ in a 
positive and negative direction with arg(£) = — tt and arg (£) = ir respectively. 
Then 



v x (x) 



2an sin (an) 



c r a dr 



< a < 1 



/o \n + e a7ri r a \' 
The asymptotics in (|4.12l) can be obtained rewriting (|4. 25[) as 



and 



v(x) = (l + a) r 



v(x) = (1 - a) + 



1 



C a + nJ ( 



C a + « J C 



-d(. 
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Differentiating these formuias we obtain 

and using the change of variables C,x = ( in both integrals it follows that 
v x (x) = - — -^ [ -d( , v x (x) = —- [ (—^ -)e c dC- 

m x L J 1 C + KX 7Tl X J-y \C + KX J 

Taking the limit x — > + in the first formula and x — > oo in the second and 
using the fact that 

m ] 1 C r (a) m J 1 T (-a) 

we obtain (gUg) , Using then (gSJ) we obtain (|4~Hi|) . (ETH)) . 
Formula (|4.15[) follows by integrating by residues. ■ 

We can reformulate the results in Theorem [T] in terms of the original func- 
tions (H (X) ,U (X) , V {X)) . In the following we choose k = 1 as a suitable 
normalization. 



Theorem 5 -For any < a < 2 £/iere exists a solution of j4-4\ ): i4-5\ ) satisfying 
%(— oo) = "H(oo) = U(—oo) — U(oo) = 0. It is given by 



"•W-'-5i/ T (?T?*)7« 



(4.26) 



■H a (X) = -^l , U a {X) = -J e^^ldZ (4.27) 
lim V a (X) = 1 + a , lim V a (X) = 1 - a . (4.28) 

X— > — oo X— >-oo 

// < a < 1 we have < 0, H a {X) > 0, U a (X) > for any lei. 

M^e Ziawe i/ie asymptotics 

dV a (X) 2a_ eaX as x ^_^ ae(Q,2) (4.29) 



T (a) 

dV a (X) 2a _ aX 



dX T (-a) 
2a 

2n 



e 



as X^- oo , aG (0, 1)U(1,2) 



H a (X) ~ ^^e 11 * as X -> -oo , a € (0,2) (4.30) 
T(a) 



H a (X) ~ --p-e- aX as X^oo , o G (0, 1) U (1,2) 
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U * PO ~ w w I TT e as * "> -oo . o e (0, 2) (4.31) 
1 (a) (a + 1) 

^ re- aX as I^co , a e (0, 1) 

r(-o) (1 - a) 

J/a = l: 

V„ (X) = 2e- eX , H a (X) = 2e x e- eX , W Q (X) = 2e~ x (l - (l + Ke x ) e" 6 *) 

(4.32) 



4.4 Solution of fl4~6l) 

Equation (|4. 6[) that approximates the behaviour of the solutions of (|3.6[) - (|3.8p 
if XE fa 1 and [/ is of order one, can be solved explicitly since 

E = X (log (U)-U)-V + log (V) (4.33) 

is preserved along trajectories. The solutions of (|4.6|) will be used to match the 
solutions of (|4.4[) . (|4.5|) obtained in the previous Section for 1 <C U -C ^. For 
this range of values V is almost constant if A — >• 0. Moreover, since we are away 
from the ODE regime described in Section I3T41 we can assume that XE in (|3.25[) 
is of order one and therefore E is of order one, whence also \V — 1| is of order 
one. 

Suppose that we consider a matching region (to be made precise later) where 
E is of order one, V is almost constant and satisfies V < 1. Then the first 
equation in (|4.6I) implies that U decreases exponentially. More preciesely, let us 
assume that U = 1 for X = X~. Then, the asymptotics of U, V are 

V~(l-o _ ) , U ~ exp (-a- (X - X~)) (4.34) 

for some a~ <E (0,1). The asymptotics (|4.34j) will be shown to match with 
suitable solutions among the ones described in Theorem [T] Moreover, notice 
that (|4.6|) implies that these asymptotics are valid as long as U remains larger 
than some small number eo that we can assume to be of order one, although 
small. The exponential decay of U in (|4.34l) implies that the transition between 
U = 1 and U — eo takes place on a length scale (X — X~) of order one. 

We now describe the solution when U is less than eo- For this purpose let 
us dehne by a_ and a+ the two roots of the equation E — —V + log (V). Then, 
in the limit A — > + , the curves given by (|4.33p in the plane (U, V) are approxi- 
mately the two horizontal lines {V = 1 — a_} and {V = 1 + a + } connected by 
a curve that is approximately vertical. Along this connecting curve V changes 
an amount of order one and Alog(J7) continues to be of order one since E is 
preserved. When U is less than eo we can approximate the second equation in 
flU} as ^ = X (1 + O (e )) V. Then V (X) increases exponentially. The first 
equation in (|4.6p indicates that for the solutions of (|4.6[) under consideration, 
U remains small while V changes from (1 — <z_) to (1 + a + ) . Then we can use 
the approximation 

V = (l- a-) exp (A (1 + O (s )) (X — X~ — i^ rans )) (4.35) 
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where (^ rans is the range of X that it takes for U to decrease from U = 1 to 
U = £o- We recall that i^ rans is of order one. 

Notice that the structure of the curve (|4.33[) implies that eventually U be- 
comes again of order Eq with V close to (1 + a + ) . It then follows from (|4.6[) that 
for small A, U reaches again the value 1 at X + after an additional length if rans 
of order one. Notice that, using (|4.35[) . we have 

V = (1 - a_) exp (A (1 + O (e )) (* + - A"" - 4n s - ^„.)) (4-36) 

and U has the behaviour, for So < U j, 

U ~ exp (a+ (X - X + )) . (4.37) 

We need to estimate the relation between a_ and a + . From (|4.33p we deduce 
that to leading order 

log (1 - a-) + a_ = log (1 + a+) - a+ . (4.38) 

On the other hand (I4.36|) allows us to obtain an approximation for the tran- 
sition length (X + — X~) . Indeed, we have 

(1 + a+) = (1 - a-) exp (A (1 + O (e )) {X+ - X~ - l+ ans - £^ rans )) 

whence, since if rans , P-trans are °^ or der one and eo can be made arbitrarily 
small 

( ^-X-)4log(l±|±) (4.39) 

as A^ 0+. 

We have described the transition of the solution (U, V, H) of (|3.6p - (|3.8[) for 
small A and [/<!. The following result will play a crucial role in describing the 
solutions, because it will show that the amplitude of the oscillations increases 
with X. 

Lemma 6 Suppose that a_ > 0, a + > satisfy \4-38\ ). Then a + > a_. 

Proof. This follows from the inequality log (1 + x) — x > log (1 — x) + x that 
is just a consequence of log (1 + x) — log (1 — x) = J Q X jt^) > ■ 

4.5 Coupling of the regimes 

4.5.1 Description of the iterative procedure 

We now describe the behaviour of the solutions of (|3.6I) - (|3.8[) for A — > and 
XE m 1, with E as in (|3.25[) . It turns out that these solutions can be described 
for this range of values by a sequence of intervals where the solutions can be 
approximated alternately by solutions of (I4.4[) . (14.51) or by solutions of (|4.6[) with 
both types of regions connected with suitable matching regions. 
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More precisely, suppose that is a point where U (X+) = 1 , V (X+) = 
1 + a n . Notice that at such a point we can expect to be able to use the approx- 
imation (|4.6[) . This will be seen matching the ODE-Integrodifferential equation 
regime with the ODE regime described in Section 13.31 We will assume for the 
moment that this approximation is valid. Then H is also of order one and we 
can approximate it to the leading order as H — UV (cf. Section 14.21) . We 
can also use the approximation (|4.37[) with a + = a n and approximate V by a 
constant. Then we obtain for \X — X+\ of order one that 

U ~ exp (a„ (X - X+)) , V~l + a n , H ~ (1 + a„) exp (a n (X - . 

(4.40) 

We can now match the asymptotics (|4.40[) with those obtained for the solu- 
tions of (|4.4p . (I4.5[) . We use (|4.1[) to obtain the matching condition 

U(X)~\exp(a n (X-X+)) , Vp0~(l + O, 

H (X) ~ A (1 + a n ) exp (a n (X - X+) ) . 

This behaviour must be matched with the one for the solutions of (I4.4[) . (14.51) 
obtained in Theorcm[5] More precisely we will match these behaviours with the 
ones of 

<f„ (X) = (U an (X - X° n ) , V a „ (X - X°) , H an (X - X°)) 

for a suitable choice of a n and X®. The asymptotics of V (X) combined with 
the one in (I4.28[) yield a = a n , since in the matching region we expect to have 
(X - X°) -> -co. On the other hand g3DJ|, (|43T|) give the choice 

X° n = Xt + - log L 2a " \) • (4.42) 
a„ \(1 + a n ) V (a„) XJ 

With these choices of a n and X° we obtain, using (|4.28|) . (|4.30|) . (|4.31|) and 
(|4.41|) , a matching to the first order between (|4.40[) and $ n (X) in the common 
region of validity where (X - X+) > 1 and (X - X°) < 1. 

We can then use the approximation <£„ (X) to describe the solutions of 
(|3.6[) - (|3.8p for small A and (X — X°) of order one. This approximation breaks 
down for (X — X^j sufficiently large. Indeed, (|4.30l) and (|4.31l) imply that 
U a (X — X°) and % a (X — X°) converge exponentially to zero as (X — X®) — > 
oo. However, the approximation (|4.4|) . (|4.5[) is only valid if \(U (X) , H (X))| > A. 
In the region where \{U (X) , H (X))\ becomes of order A we must use again 
the approximation (|4.6|) . In particular this region can be described using the 
analysis in Section 14.41 The asymptotics of the solutions for U of order one is 
as in (|4.34l) for suitable choices of a_, X~ . More precisely, since this matching 
is made for (X — X°) > 1 we can use the asymptotics (|4.28[) . (|4.30p . (|4.31l) for 
I^oo. We then obtain, using also the rescaling (|4.1[) . the matching conditions 

t^W~- r , wT ^-e— (*-^), F(X)~(l-a n ) 

1 {—On) (1 - a„) A 

(X) ~ - 2a- e -^(x-^) for (x _ x o } 1; ^ x 
r(-a„)A 
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We now match these formulas with (|4.34l) (combined with the approximation 
H = UV). To this end we must choose 

a-=a n , X- = X° n + - log (-- ^ -\) . (4.43) 

a n \ (1 - a n )T (— a„) A/ 

We remark that all this analysis will be meaningful only if < a n < 1. 
Therefore T (-a n ) < 0. 

The region where U < 1, H < 1 can then be described using the analysis 
in Section l4~4l The conclusion of this analysis is that (U, V) moves close to the 
point (1, (1 + a n+ i)) for some suitable a n +\ in a characteristic length given by 
P~3U]) . More precisely if we define (cf. (0735])) 

log (1 - a n ) + a n = log (1 + a„+i) - a n+1 , a n > , a n+1 > (4.44) 

X„ ++1 = X- + Il„ g (l±^±l) (4.45, 

we can approximate {/, V by means of (|4.40l) with replaced by 
Combining (|4.42l) , (|4.43p and (|4.45l) we obtain to the leading order 




200 400 600 2 4 6 8 10 

X u 



Figure 4.1: Intermediate regime: U (left) and phase plane for U, V (right) 



4.5.2 Matching with the ODE regime 

We first remark that equation (|3.19[) agrees with the approximating equations 
(14.41) . ()4.5|) and (|4.6j) if E 3> 1, i? <C j in their respective regimes of validity. 
Indeed, we can rewrite (|3.19[) using (|3.18l) as 

% = U{V-l) + \U{U-l) , % = \{l-U)V (4.47) 
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and these equations yield the same behaviour as (|4.6[) if U is bounded. 

On the other hand, in order to approximate the solutions of (|4.4p . (14.51) for 
£> 1, £ « j we use the fact that for small a, the solutions of (|4.4I) . (|4.5p 
in Theorem [5] have the characteristic length scale - . Then, for a small, which 
corresponds to the matching region indicated above, we have a characteristic 
length very large compared with one. Then, as in (I3.14p we can approximate 
I [H] ~ XHU. Since E < \ we have U (X) < 1 and we obtain from (|4~4l 
the approximation % (X) = W (X) V (X) + (W (X)) 2 V (X) and plugging this 
formula into (|4.5|) and using the rescaling (14.1[) we obtain 

g = £/ ( U-l) + A([/) 2 U , ^-A[/U-A 2 (L/) 2 U 

that agrees with (|4.47l) in the range of values U ^> 1, [7 <C t • 

The previous computations show that the equations used in both transition 
regimes agree in the intermediate matching regime. Actually it is possible to 
check that the main characteristics of the computed solution also agree, as could 
be expected. Notice that the length scale {X^+i ~ -^n) can ^ e computed using 
the function T (E) in (ET24]) for E > 1, 

T(E)~2y/2VE. (4.48) 

In the computation of this integral, we have split the region of integration 
into the two subsets (J7_ (E, 0) , sq) and (eq 7 U + (E, 0)) where Eq > is a small 
but fixed number. It turns out that the contribution of the second integral 
is much smaller than the one due to the first integrand, exactly in the same 
manner that the contributions in (X^ +1 — X£) computed in Subsection 14.5.11 
are due mostly to the region with small values of U. Using the rescaling (|3 . 18[) 
we then obtain the approximation 

OCh - - ■ (4-49) 

On the other hand, we can compute the same length using (I4.46[) . Taking 
into account that in the matching region a n , a n+1 are small, we obtain from 
(gH) that 

2 2 

a n+1 ~ a n + - (a n ) . (4.50) 
Using (I4.46P we then obtain to leading order that 



x: +1 x: ~ - ^ - =f ■ (4 - 51) 

Next, equation (|3.25[) and the definition of a n by means of U = 1, imply 

a n = V2VeV\ (4.52) 
and plugging (|4.52l) into (|4.51l) we obtain the sought-for matching with (|4.49|) . 
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Finally we obtain the matching for the formula of the change of energy 
in each iteration. We have obtained using the ODE approximation that the 
change of energy in each cycle is given by (|3.26[) where in the matching region 
1 <C E j we can use the approximation (I3.28|) . Then 

E n+1 -E n ~ -^-E% . (4.53) 

On the other hand, using again a n = \f2\f~E^f\ and plugging it in (|4.50p we 
obtain again (|4.53|) and this yields the desired matching. 



5 Shooting argument 

The structure of the solutions of (|3.6[) - p.8[) described in Section POl bv means of 
alternate approximate solutions of (|4.4[) - (14.51) and (|4.6[) yields several sequences 
of numbers {o„}, , {X~} , {X®} that measure the amplitude of the 

oscillations, the positions of the points where U = 1, V > 1, U = 1, V < 1 
and the position of the peaks where U and H are of order \ respectively. It is 
relevant to notice that the sequence {a n } is increasing due to Lemma [5] The 
values of a n approach zero in the matching region (see the approximation f|4. 52[) 
for 1 <C E <C j), but they become of order one for E « j- with increases in 
each iteration of order one. Eventually, this sequence of points becomes larger 
than or equal to one. 

A detailed study of the solutions of (|3.6l) - (|3.8p shows that there are several 
possibilities. If a n becomes strictly larger than one, Theorem [5] shows that 
during the peak regime V becomes strictly negative. In particular V reaches 
the value V = at some finite X*. If V remains positive for larger values of 
X, something that could happen if a n = 1, it would be possible to approximate 
(|3.6l) - (|3.8p by means of (|4.6p . Therefore V would increase, and the value of a n 
would be increased to a new value a n +i > 1 at the next iteration. The only 
possibility of not having such a behaviour would be with V — > as X —> oo. 
Such approximation to zero cannot be described by the approximate equation 
(|4.6p . because in such a regime the integral terms in (|3.6[l - (|3.8[) must be taken 
into account in full detail. 

We now recall that for any K £ C there exists a unique solution of (|3.6p ~ 
(|3.8p satisfying p. lip . The discussion in Section 13.21 shows that all the so- 
lutions of (|3.6p ~ p.8p with such a behaviour can be obtained, up to rescaling, 

by choosing K in the interval l,exp • A continuity argument tak- 

ing into account the behaviour of the trajectories shows that there exists a 
value of K — K* in this interval for which the corresponding solution of (|3 . 6|) - 
(|3.8p satisfies limx^oo V (X) = 0. This also implies limx-j-oo H (X) — and 
limx^oo U (X) = using (|3.6p . (|3.7p . More precisely, the inversion of (|3.6p 
allows us to write H (X) in terms of UV. Then limx^oo H (X) = and (|3.7p 
yields lim^oc U (X) = 0. 
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6 Asymptotics as X — > oo 

We finally describe the asymptotics of the solution of (|3.6[) - (13.81) satisfying 
(|3.1ip with K = described in the previous Section. We have seen that 
lim^^oo U (X) = limx^oo V (X) = limx-s-oo H (X) = . To leading order the 
behaviour of this trajectory is described by the solution of (|4.4I) . (|4.5I) given in 
Theorem [5] with a = 1. The asymptotics 

iT(X)~^ 1 (6.1) 

is valid as long as H ^> 1, (X — X®) ^> 1, where, by assumption, to the leading 
order a n = 1. It is convenient to reformulate (16.11) in the original variables (cf. 
dOt . (1331) ) 

2 / x \ / / x 



(ar) ~ — ( — ) exp ( — (__) ) (6.2) 



that is valid for f- > 1, ft > 1. 

We can match (|6.2[) with the following exponential asymptotics for the so- 
lutions of (|2.7[) as x — > oo: 

2 

/i (x) ~ t- (kx) exp (-ki) (6-3) 
A 

The asymptotics (|6.2j) . (|6.3j) match for small A if K = — . This gives the 
desired exponential decay for the obtained solution. 



7 Self-consistency of the approximations of / [H] 

We now show that the solution we obtained is self-consistent with the approxi- 
mations made for the integral operator / [H] (X) defined in (|3.8[) . We have made 
two main approximations. The first one is (|3.12p . that approximates I [H] by an 
operator with a simple dependence on A. The second approximation is f|3 . 14[) and 
it allows us to replace the integral operator by a much simpler local operator. 

Concerning the validity of (|3.12[) we remark that its precise meaning is that 
the error made in the approximation is smaller than the right-hand side. In 
order to check its validity we distinguish between the regions where H is of 
order 4 and the regions where H is smaller than that quantity. 

Notice that for any region we can expect the errors made in the approxima- 
tion to be of order 

dYeV- x) (X-Y)H (Y) / dZH (Z) , 

-oo J\og(e x — e Y ) 

R 2 - A 2 [ X dYe {Y - x) H (Y) C dZ (X - Z) H (Z) 

J~ oo J\og(e x —e Y ) 
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and we can expect these two terms to be very small compared with the right- 
hand side of (|3.12p . Indeed, in the case of Ri, the contribution due to the region 
A (X — Y) < Sq is small compared with the term in (13. 12)) if Eq is small. On the 
other hand, if A (X — Y) > So we can use the smallness of the exponential factor 
e (y—x) to obtain estimates for the corresponding terms in R\ as CX 2 e~~*~ 
and since H-H]^ < C/X it then follows that this contribution is exponentially 
small. Concerning R2 we can argue similarly for A (X — Z) smaller than £0. If 
A (X — Z) > £ , since Z > X + log (l - e Y ~ x ) it follows that the size of the 
region of integration can be bounded by Ce~~^ and therefore, it gives also a 
very small contribution. 

Concerning the approximations of / [H] by means of local terms that have 
been made in the derivation of the ODE approximations (j3. 19[) or (14.61) two 
main ingredients are required. In the derivation of (|3.19[) we have used just the 
fact that H (Y) has significant changes over distances much longer than one. 
Since the order of magnitude of H is roughly the same for the range of values 
described by means of p.!9[) no special care is required to estimate the values 
of H (Y) with X — Y 1. On the other hand, in the case of the approximation 
(|4.6p we use the fact that I [H] is a corrective term that is completely ignored 
in (|4.6[) . Notice that H (Y) takes values for Y < X much larger than the ones 
of H (X) in the region where the approximation (|4.6p is used. However, due to 
the exponential factor e^ Y ~ x ^ as well as the exponential decay of H (Y) in the 
region described by the integro-differential equation (cf. Subsection I4.5.1[) the 
corresponding contribution of such values of H (Y) is negligible. 



A Uniform estimates for the change of energy 

During the part of the dynamics dominated by a perturbation of a conserva- 
tive ODE it has been assumed that the energy E does not change significantly 
during each cycle of length T (E) . This is the justification for the adiabatic 
approximation and we will check now that, indeed, the adiabatic approximation 
yields such smallness for the variation of the energy during the trajectory and 
thus the assumption is self-consistent. 

Let us denote by ^0 € t a value where U (£0) = 1, V (£0) > !• Our goal is to 
estimate this quantity for £ < £ < £ + T (E (£ )) , where £ + T (E (£ )) is the 
next value of £ where U (£) = 1 and V (£) > 1 again. Due to the invariance of 
(|3.19p we can assume without loss of generality that £0 = 0. Due to (|3.23p and 
assuming the adiabatic approximation we obtain the following approximation 
for the change of the energy along a trajectory: 



E (Q - E (0) 



(U (s) - if + (oj (s)) 2 (1 - U (s)) ds = a(0- (A.l) 



Since 



^Jf- = (U (0 - 1) \(U (0 - 1) - (wo (O f 
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and given the form of the curves (Uo, uiq) satisfying p.22[) we obtain the existence 
of numbers k = 1, 2,3 satisfying = £ < 6. < 6 < & < 6)+r (-B (Co)) = & 
as well as 

^1 = 0, fe = 0,1,2,3,4. (A.2) 



Moreover [U (Co) - 1) = (U (&) - 1) = (Co (60 - 1) - 0, v/t/ (6)-l 
w (60 = -^0 (6) and 

" rr! °<0 if £e (6,60 u (6,60 , ^#>o if £e(6,6)u(6,6) 



(A.3) 

The value of the total change of the energy during each cycle, namely a (6) , 
has been computed in (I3.26[) . (I3.27[) and approximated asymptotically as E — > oo 
(cf. p~28|) ). Notice that a (6) =$(£'), assuming that £ (6) = 

Due to (|A.2j) . (|A.3j) we need to compute cr (6) , cr (6) , cr (6) in order to 
estimate the range of variation of a. We then have 



>(6) , ^ (6)} < <7 (0 < m ax {a (6) , a (£4)} 



(A.4) 



where we use that a (6) > a (£1) . In order to estimate these quantities we define 
an auxiliary quantity f2o via E = JIq /2 that measures the maximum value of Wo 
for a given value of E. It will be more convenient to compute cr (£) using wo as 
independent variable. To this end we define U— (loq; Qq) < 1 < U+ (ujq; fio) by 
means of the roots of the equation 



(tf-l)-]og(E0 + y = ^ 



(A.5) 



Taking into account p. 211) we can use wo as variable of integration instead of 
s in order to compute cr (£) for £ <E (0, £3) . Then, with u>k := luq (6) we obtain 



c (60 = crfc (^0) 



(^+(wo;n )-l)-(wo) dw o , fc = l,2,3, (A.6) 



where the integration is made along the curve defined by (|A.5|) . Notice that we 
use w (£3) = -^o- 

In order to approximate cr (61) , cr (6) , cr (£3) for small fio we use Taylor in 
(|A.5|) to obtain the approximation 



(D±(w ;no)-l) a +wg = ng 



(A.7) 



that combined with (|A.6[) yields 



O" (6c) = CTfe 



\M ~ w o - ("b) S 



dwo as J7o — > . 



On the other hand, we can approximate wi, CJ2 by ^/ — u;^ = {usk) 2 for 
fc = 1, 2 whence cj^ ~ ±^o as Qq — > for k = 1, 2 and 



(A.8) 
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All these quantities must be compared with a (£4) that can be computed 
using (|3.27[) . The approximation (IA.7|) then yields 

cr (£4) ~ ttQI as fi -> . (A.9) 

Combining (|A.8[) . (|A.9|> . as well as (|A.4[) it follows that 

o (fig) < ff (0 < wfio 1 f e (°> &) and ^o^0+. 

We can obtain a similar estimate for large values of Qo- Using (|A.5[) we 
obtain 

U+ (w ;O ) - 1 ~ - (fig -o>o) as ft 00 

that is valid as long as (U — 1) is large. 

Since the region where (U — 1) is of order one gives a small contribution to 
the integrals, we obtain the approximations 

/3 

a = = -<r(&) , <r(&)=o(fi§) as^o-^oo. 

Combining this with (|3.28[) we obtain that 

-CQl < a (£) < , £ G (0, £4) and ft -> 00 

for some C > independent of r2 - It then follows that the change of energy 
during each cycle can be estimated by the final change. Therefore, this justifies 
the adiabatic approximation. 

Acknowledgment: BN and JJLV gratefully acknowledge the warm atmo- 
sphere at the Isaac Newton Institute for Mathematical Sciences where part of 
this work was done during the program on PDE in Kinetic Theories. This work 
was also supported by the EPSRC Science and Innovation award to the Ox- 
ford Centre for Nonlinear PDE (EP/E035027/1) and through the DGES Grant 
MTM2007-61755. 



References 

[1] M. Abramowitz and I. A. Stegun, Handbook of mathematical functions with 
formulas, graphs and mathematical tables, National Bureau of Standards 
Applied Mathematics Series, 55 (1964) 

[2] J. Bcrtoin, Eternal solutions to Smoluchowski's coagulation equation with 
additive kernel and their probabilistic interpretation, Ann. Appl. Probab. 
12 , (2002), 547-564. 

[3] J. Cahizo and S. Mischler, Regularity, asymptotic behavior and partial 
uniqueness for Smoluchowski's coagulation equation, to appear in Rev. Mat. 
Iberoamericana, (2010) 



26 



[4] R. L. Drake, A general mathematical survey of the coagulation equation, 
In: Topics in Current Aerosol Research (part 2). International Reviews in 
Aerosol Physics and Chemistry, Oxford. Pergamon Press (1972), 203-376. 

[5] M. Escobedo, S. Mischler and M. Rodriguez Ricard, On self-similarity and 
stationary problems for fragmentation and coagulation models, Ann. Inst. 
H. Poincare Anal. Non Lineaire 22, (2005), 99-125. 

[6] M. Escobedo and S. Mischler, Dust and self-similarity for the Smolu- 
chowski coagulation equation, Ann. Inst. H. Poincare Anal. Non Lineaire, 
23, (2006), 331-362. 

[7] F. Filbet & P. Laurengot, Numerical simulation of the Smoluchowski equa- 
tion, SIAM J. Sc. Computing, 25, (2004), 2004-2028. 

[8] N. Fournier & P. Laurengot, Existence of self-similar solutions to Smolu- 
chowski's coagulation equation, Comm. Math. Phys. 256, (2005), 589-609. 

[9] N. Fournier & P. Laurengot, Local properties of self-similar solutions to 
Smoluchowski' s coagulation equations with sum kernel, Proc. Royal Soc. 
Edinb. 136 A, (2005), 485-508 

[10] P. Laurengot and S. Mischler, On coalescence equations and related mod- 
els, Modeling and computational methods for kinetic equations, Eds. P. 
Degond, L. Pareschi, G. Russo, Series Modeling and Simulation in Science, 
Engineering and Technology, Birkhauser, (2004), 321-356. 

[11] M. H. Lee, A survey of numerical solutions to the coagulation equation, J. 
Phys. A 34 , (2003), 10219-10241. 

[12] F. Leyvraz, Scaling theory and exactly solvable models in the kinetics of 
irreversible aggregation, Phys. Reports, 383 2/3, (2003), 95-212. 

[13] G. Menon and R. L. Pcgo, Approach to self- similarity in Smoluchowski's 
coagulation equations, Comm. Pure Appl. Math., 57 , (2004), 1197-1232. 

[14] B. Nicthammer and J. J. L. Velazquez, Optimal bounds for self-similar 
solutions to coagulation equations with product kernel, Comm. PDE, to 
appear (2011). 

[15] B. Niethammer and J. J. L. Velazquez, Rigorous construction of a self- 
similar solutions to a coagulation equation with multiplicative kernel, work 
in preparation. 

[16] M. Smoluchowski, Drei Vortrage uber Diffusion, Brownsche Molekularbe- 
wegung und Koagulation von Kolloidteilchen, Physik. Zeitschrift 17, (1916), 
557-599. 

[17] P. G. J. van Dongen and M. H. Ernst, Scaling solutions of Smoluchowski's 
coagulation equation, J. Stat. Phys., 50, (1988), 295-329. 



27 



